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Abstract. We investigate certain configurations of Luneburg lenses that form light 
' ^ . propagating and guiding networks. We study single Luneburg lens dynamics and apply 

^-j" the single lens ray tracing solution to various arrangements of multiple lenses. The 

*-£h . wave propagating features of the Luneburg lens networks are also verified through 

direct numerical solutions of Maxwell's equations. We find that Luneburg lenses may 
form efficient waveguides for light propagation and guiding. The additional presence 
of nonlincarity improves the focusing characteristics of the networks. 



PACS numbers: 41.20.Jb,42.15.Dp 
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Hamiltonian ray-tracing, waveguides, Kerr effect. 



2 



1. Introduction 

Gradient Index (GRIN) metamaterials are formed through spatial variation of the index 
of refraction and lead to enhanced light manipulation in a variety of circumstances. 
These metamaterials provide natural means for constructing various types of waveguides 
and other optical configurations that guide and focus light in specific desired paths. 
Different configurations have been tested experimentally while the typical theoretical 
approach uses Transformation Optics (TO) methods to cast the original inhomogeneous 
index problem to an equivalent one in a deformed space [U [2J El IU [6]. While this 
approach is mathematically elegant, it occasionally hides the intuition obtained through 
more direct means. Furthermore, a general, continuous GRIN waveguide may be hard 
to analyze in more elemental units and relate its global features to these units. In the 
present work, we adopt precisely this latter avenue, viz. attempt to construct waveguide 
structures that are seen as lattices, or networks, of units with specific features. This 
is a "metamaterials approach", where specific properties of the "atomistic" units are 
inherited as well as expanded in the network. 

The "atomic" unit of the network is a Luneburg lens (LL); the later is a spherical 
construction where the index or refraction varies from one, in its outer boundary, to \/2 
in the center with a specific functional dependence on the lens radius [HE]. Its basic 
property, in the geometrical optics limit, is to focus parallel rays on the spherical surface 
on the opposite side of the lens [H El El [91 HUl HH [12l H3l [H] . This feature makes LL's 
quite interesting for applications since the focal surface is predefined for parallel rays 
of any initial angle. While the rays traverse the lens, they suffer variable deflections 
depending on their distance from the optical axis ray, leading to a point image on the 
lens surface. This property of the spherical LL is also shared by its cylindrical equivalent 
formed by long dielectric cylinders while the light wave-vector impinges perpendicularly 
to the cylindrical axis. This geometry turns the problem into a two dimensional one, 
constructed for any plane that cuts the LL cylinder perpendicular to its axis. The 
work in this article focuses on exactly this type of cylindrical LLs and, as a result, our 
approach is strictly two dimensional [7J El [10] . 

In the present work we analyze the ray trajectories in a single LL and subsequently 
we use this analysis in waveguides build exclusively through LLs. Specifically, in 
the next section we present analytical approaches for the trajectories traversing a LL 
[H El El [TJ El Uni HS1 [LEI mj [HI [19] and derive explicitly a solution expressing the ordinate 
of the trajectory as a function of the abscissa in Cartesian coordinates. In section 3 
we introduce the main point of the present work, viz. the formation of Luneburg Lens 
Waveguides (LLW) through the arrangement of multiple LLs [13] in geometrically linear 
or bent configurations. The LLWs are investigated through dynamical maps stemming 
from the analytical solution for the single LL. In section 4 we investigate the same LLWs 
but using a direct solution of Maxwell's equations and find compatibility with the results 
of ray tracing approach of the previous section. Additionally, we also use nonlinearity 
in the index of refraction and find better focusing characteristics [121 120] • Finally, in 
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section 5 we conclude. 



2. Single Luneburg lens 

2.1. Ray tracing solution 

The general problem we have at hand is light propagation in an inhomogeneous isotropic 
medium with an index of refraction n(r), where r is the position vector. Although this 
problem has been tackled though various approaches, for the purposes of our problem 
we need the exact solution for propagation through a single LL, which is itself an 
inhomogeneous medium. What we need is a ray tracing solution for a single LL that, 
in turn, will be used for providing the basic element in a map approach for the ray 
dynamics in LL networks. 

Given the recent interest in LLs we present below three ways for obtaining the ray 
tracing solution for the single LL; these explicit solutions from the different approaches 
may be useful in other LL studies as well. The first two methods (section 2.2 and 
Appendix A.l respectively) are based on the Fermat principle for optical path or travel 
time minimization [Q [7J [151 EH EH EB] while the last one (Appendix A. 2) uses a 
geometrical optics approach to Helmholtz wave equation [U [7J, [T7J, [T8j UHl [23] . In all cases 
we focus on the specific, cylindrically symmetric, LL index of refraction n(f) [7l[8l [T0l[T9] 

<r) = ^2 -(I)' (1) 

where R is the radius of the lens while r is the radial distance from the center in the 
interior of the lens, i.e. r ^ R. The LL is supposed to be embedded in a medium with 
index equal to one, leading to a continuous change in the index variation as a ray enters 
the lens. 



2.2. Quasi 2D ray solution 

The optical path length S of a ray from a point A to B is given by [U \T5\ dSl ETJ Q2] 
S = f nds (2) 



J A 

In polar coordinates the arc length is ds = yl 1 + r 2 <j) 2 dr, with = dtp/dr where 
the coordinate r is considered a "generalized time". As a result Fermat's variational 
integral of Eq. ([2]) , becomes 



S= I n(r)y 1 + r 2 (f) 2 dr (3) 

J A 

and leads to the Optical Lagrangian [TJ [151 HH] 



£(<f),<j),r) = ra(rWl + r 2 2 (4) 



4 



The shortest optical path is obtained through minimization of the integral of Eq. and 
may be found through the solution of the Euler-Lagrange equation for the Lagrangian 
of Eq. (ID, viz. 

d dC dC 

Since the Lagrangian (jlj) is cyclic in <j), we have d£/d<fr = and thus dC/dcf) = C where 
C is a constant. We thus obtain for the Lagrangian of Eq. (jlj) [TJ \TE\ [TB] 

n(r)r 2 



J 1 + r 2 <p 2 



C (6) 



This is a nonlinear differential equation describing the trajectory r(0) of a ray in the 
interior of LL. Replacing the term = d<fi/dr and solving for dip, we obtain a first 
integral of motion, i.e. 

fdr (7) 



r\fn 2 r 2 — C 2 

The expression of Eq. ([7]) holds for arbitrary index functions n(r). Using the specific 
LL refractive index function of Eq. (CD) we evaluate the integral (I7|) and obtain the ray 
tracing equation for a single LL as follows 

C'R 



^ 1 - y/T^C# sin(2(0 + 0)) 



(8) 



where C and (3 are constants. This analytical expression may be cast in a direct 
Cartesian form for the (x, y) coordinates of the ray; after some algebra we obtain 



(1 — T sin(2/3))x 2 + (1 + T sin(2/3))?/ 2 

-2T cos(2(3)xy + (T 2 - l)i? 2 = (9) 

where T and /3 are constants. We note that Eq. Q is an equation for an ellipse. This 
result agrees with the Luneburg theory and states that inside a LL light follows elliptic 
orbits. 

The constants T and /3 of Eq. are determined through the ray boundary (or 
"initial" conditions) and depend on the initial propagation angle 9 of a ray that enters 
the lens at the point (xo,yo) lying on the circle at the lens radius R [TS]. The entry 
point of the ray is at (x,y) = i?(cos(6> + ir), sin(6 ) + %)) = —R(cos9,sin9). Substituting 
these expressions in Eq. ((9]) we obtain after some algebra 

T = sin(2/3 + 29) (10) 

In order to determine both constants T, /3 we need an additional relation connecting 
them. Taking the derivative of the Eq.()9]) respect to x and regarding the fact that 
dy/dx = tan(6*), where 9 the initial propagation angle. In addition, using the labels 
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(x ,y ) for the initial ray point on the LL surface, we set x = x and y = y in Eq.flHj) 
we solve for T and obtain 

T= x + y tan(fl) 

tan(0) [x cos(2/3) - y sin(2/3)] + [x sin(2/3) + y cos(2/3)] 1 J 

The Eqns. fflUj) and ffTTj) consist of an algebraic nonlinear system expressing the 
constants T and (3 as a function of the initial ray entry point in the LL at (xo, yo) with 
angle 9. Combining the Eqns. ( flOj) . (fTTj) and after simplifications we obtain 

/3 = i(tan- 1 (x o /y o )-0) (12) 
therefore, according Eq. ffTO]) 

T = sin(tan _1 (xo/2/o) + 0) (13) 
Substituting now the Eqs. (ll2|ll3l) to Eq.(j5j) and solving for y, we obtain the equation 

(2x y + R 2 sin(29)) 
V[X) 2x 2 + (l + cos(26))R 2X 

V2Ry cos(#) ^(1 + cos(29))R 2 + 2x\ - 2x 2 
+ 2x 2 + (1 + cos(20))f2 2 

x sin(0) v/(l + cos(20))i? 2 + 2x^ - 2x 2 



(14) 



2x 2 + (1 + cos(29))R 2 

where Xo>2/o gi ye the initial ray entry point and 9 the initial propagation angle. 

The Eq. (|14|) describes the complete solution of the ray trajectory through a LL; 
in the simplest case where all rays are parallel to the x axis, the initial angle 9 = and 
Eq. ( fl4|) simplifies to 



= -2^( x ox + rJr 2 + x 2 - x 2 ) (15) 
Xq + R z v 

We note that in order to determine the exit angle 9', i.e. the angle with which each 
ray exits the lens, we may take the arc tangent of the derivative of Eq. (114j) wrt x at 
the focal point on the surface of lens at x = Rcos(9). The solution of Eq. (I14p may be 
used to study several configurations of LLs; however special provisions are necessary in 
cases where the rays need to turn backwards. It is thus more practical to use parametric 
solutions where the ray coordinates x, y are both dependent variables. This approach is 
explained in the Appendix where the parametric solution of Eq. (IB. 161) is derived. 

3. Luneburg Lens Waveguides 



The analytical solutions for the single LL trajectory, viz. Eqs. (fT4|) or (1B.16j) may be used 



in order to study analytically the ray transfer through various configurations of LLs that 
form waveguides [13]. Using the initial entry point on the LL (xq, yo) as well as the initial 
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Figure 1. The red dotted lines denote the arrangement of lenses. The blue lines show 
the ray tracing. Light is guided by Luneburg lenses across the linear network, (i) 
Arrangement with 7 lenses (left figure). We can see that all rays are focus on the last 
lens, (ii) Arrangement with 6 lenses (right figure). All rays exit in the same mode as 
originally entered, i.e. parallel to the waveguide axis. 

ray angle 9 we obtain through Eq. ( 1T4|) the exit point (x, y) and the associated exit 
angle 9'. We may thus form a mapping from (x , yo, 9) to (x, y, 9'); further propagation 
in the surrounding medium is rectilinear while the entry to the next LL is governed by 
a new initial entry point with angle equal to the previous exit angle. The resulting ray 
may be traced easily through the map. 

A geometrically linear arrangement of touching LLs on a straight line is shown in in 
Fig. [T| this configuration forms a waveguide that channels the light through. Depending 
on the number of lenses even or odd number) the rays focus in the last LL surface or 
exit as they entered respectively. In both cases we sent a beam of rays parallel to the 
axis of symmetry of the network. 

In Fig. (T5]) we use the ray tracing equations, Eqs.( IT4|) or ( IB. 161) to study the 
propagation along curved arrangements of LL's. We observe the efficient channeling of 
the rays through the network that leads to a complete ray banding at a right angle. 
Some rays escape in the sharp bend, but generally the guiding is very efficient for such 
a drastic change of direction. 

In Fig. ([3]) we form a full circle bend waveguide through a sequence of LLs and 
follow the light propagation in the geometric optics limit. We find that light may 
propagate efficiently through a loop, signifying that arbitrary waveguide formation and 
guiding is possible. 

The LL network cases presented (linear, right angle, curved) signify that LLs may 
be used as efficient waveguides. Their advantage over usual dielectric guides is that light 
bending occurs naturally through the LL properties while the outgoing light may be also 
focused, if so desired. In bends there are naturally some losses that, in the geometric 
optics limit, may be estimated by comparing the number of the incoming to the outgoing 
rays, viz. iV/jv) versus Nout respectively. In the linear arrangement of waveguides, as 
in Fig. [U the performance is perfect since N IN = N OUT . In the bend cases, such as 
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Figure 2. A right angle waveguide is formed through a network of eleven Luneburg 
lenses. The ray-tracing is performed analytically using Eq. (|B.16[) . We observe efficient 
guiding with a total bend of 90° degrees. 




X 



Figure 3. A full circle bend waveguide formed through seventeen Luneburg lenses. 
The ray-tracing is performed through the ray tracing methods. We see that Luneburg 
waveguides can guide light in a complete circle bend trajectory; more complex paths 
are also possible. 

in the right angle arrangement of Fig. [2] as well as in the full circle bend waveguide 
of Fig. [3j we find Nout/Nin = 13/21 = 0.6. We note that the aforementioned losses 
depend also on the ray coverage of the initial lens as well as the sharpness of the bend; 
the losses can be reduced by manipulating appropriately these two factors. 
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4. Wave propagation in Luneburg lens networks 

4-1. Propagation through LL waveguides 

In the present section we present results for the true wave propagation in networks of 
LLs. The simulations we done both with a homemade FDTD based algorithm that was 
tested against COMSOL Multiphysics® software. Both gave good agreement and, as a 
result, we show only the results obtained from the latter. In order to be more realistic, 
we used Gaussian beams instead of plane waves [T2| [20] . In the single LL case the effect 
of the gaussian beam shape is to shift the focus outside the LL surface, an effect that 
may be compensated though nonlinearity [T2"j I2U]. Specifically, we simulate all three 
arrangements shown previously, viz. Figs. (TJ [2] and [3] and we calculate numerically the 
losses of the guiding. Subsequently, we introduce nonlinearity and show that the Kerr 
effect improves guiding as well as focusing of the beam through Luneburg waveguide. 

In order to setup our arrangements we use silica glass for the Luneburg lenses with 
refractive index variation as in Eq.flTJ. In addition we use air for the outer medium 
with index of refraction n — 1 [TTJ (T2J [HJ E] with a monochromatic TM mode Gaussian 
EM wave. We perform all simulations in microwave regime with wavelength A = 1cm. 
In all simulations the radius of each LL is ten times greater that the wavelength, i.e. 
R = 10A = 10cm. Finally, in the lattice edge we apply PML boundary conditions [2~T] . 

In Figs. (|4}17j) we show the results for the full propagation of EM waves through 
LLs waveguides. We plot the steady state of the absolute value of the electric field. 

For the linear arrangement of LL waveguides we study a sequence with seven and 
six lenses, Figs. |H and |5] respectively. In the case with even number of lenses, as in Fig. 
HI the ray beam is guided and defocusses while exciting the waveguide. On the other 
hand, in the case with odd number of lenses, as in Fig. El the beam focuses in the last 
lens. The results obtained are compatible with those obtained through the ray tracing 
methods of Fig. [T] The right angle bend is shown in Fig. El while in Fig. [7] we have 
a full circle arrangement of LLs; these arrangements are the same as those analyzed in 
Fig. |2] and Fig. |3] respectively. In all cases studied, the numerical solution of Maxwell's 
equations is compatible with the findings obtained through the ray tracing map. 

4-2. Propagation losses 

In order to find how efficient is the guiding by LL waveguides we may calculate losses 
due to propagation through LL networks. Subsequently, we calculate numerically the 
absolute value of the electric field in the first and in the last lens in each arrangement. 
Specifically, we calculate the integrated field intensity at a cross section that is 
perpendicular to the propagate direction and is located in the center of each lens and 
compare outgoing to incoming waves. We find that in the linear case of Fig. 141151 
that Equt/Ein — 0.98. As a result, linear guiding is almost perfect. For the right 
angle bend waveguide of the Fig. El we find that Eout/Ein = 0.73 while in the circular 
arrangement of the Fig. ^Equt/Ejn = 0.70. We conclude that guiding is quite efficient 
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Figure 4. Linear Luneburg waveguide arrangement with seven lenses. Beam guiding 
and focusing on the surface of the last lens is observed. The result is in agreement 
with Fig. QJi). Furthermore, no losses are observed. 




12 Tl 2t'5»..t 



Figure 5. Linear waveguide arrangement formed with six Luneburg lenses. The EM 
wave is guided through the linear network of lenses and the wave defocusses in the last 
lens. The result is in agreement with Fig. [TJii) . Additionally, the guiding is almost 
perfect since the losses are practically zero. 



for such sharp bends. We note that losses can be reduced further through using LLs 
with different radii as well as more efficient geometries that have less sharp curvatures. 

4-3. Nonlinear Kerr effect in LL waveguides 

As already mentioned, if the propagating wave is a Gaussian beam the focus point 
shifts outside of a single LL lens. To restore focusing one may use a nonlinear medium 
[121 ED] - We use this approach in the LL waveguides as well. In the Fig. [8] we propagate 
a Gaussian beam through a linear Luneburg waveguide which is formed by three lenses 
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Figure 6. A right angle bend waveguide is formed by eleven Luneburg lenses. The 
EM wave is guided and bends at right angle through a network of lenses. The losses 
in this case are approximately 27% of the initial amplitude, i.e Equt/ Ein = 0.73. 




Figure 7. A full circle bend waveguide is formed by seventeen Luneburg lenses. 
The EM wave is guided through a circular orbit. The losses in this arrangement are 
approximately 30% since Equt/ Ejn = 0.70. 



in the absence of nonlinearity; we observe that the focus is shifts to the exterior of the 
last lens. 

In the presence of the Kerr effect the refractive index is a function of the electric 
field E, viz. 

n{E) =n L + X \E\ 2 (16) 

where x is the nonlinear factor with units m 2 /V 2 and ni is the linear part of the 
refractive index [El QH QS1 ED] • 

For the numerical simulations we use a nonlinearity factor that makes the nonlinear 
part of the index of refraction, viz. the term xl-^| 2 ) to be of the order of 1% of the linear 
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Figure 8. Linear waveguide arrangement is formed by three lenses. Simulate the 
propagation of the EM wave without any nonlinearity (Kerr effect). We can see that 
the focus points are shifted out of the surface of lenses. 




Figure 9. Linear waveguide arrangement is formed by three lenses. We simulate the 
propagation of the EM wave including the Kerr effect. We can see that the focus points 
are located again on the surface of lenses. 



part riL. In Fig. [8] we show the shift in focusing of a LL network without nonlinear term, 
while in Fig. [9] the same arrangement but with the addition of nonlinearity. We note 
that nonlinearity shifts focusing in the LL network to the surface of the lenses [T2| I2U] . 



5. Conclusion 



Luneburg lenses provide simple units for constructing GRIN-based metamaterials. The 
advantage of using LLs is that one may resort to an "atomistic" picture where the lens 
is a single unit that has specific properties. Subsequently, these properties may be used 
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in forming extended networks. The basic feature of LLs is their property to have a 
prescribed focal surface for parallel rays that lies on the lens surface. This property was 
used in the present study for the formation of waveguides made of Luneburg Lenses. 
We saw that these guides are quite versatile since they can channel light in an efficient 
way across different paths. 

In this work we obtained an exact quasi-two dimensional solution for light trajectory 
through a single cylindrical LL that can be used to study propagation in arbitrary 
lens arrangements. This solution is parametrized on the initial entry point of light on 
the surface of a LL as well as the initial ray angle and gives the output location and 
the exit angle. One may thus form a type of input-output mapping that enables the 
study of propagation in various LL configurations. While this solution is analytical, it 
fails for backward propagation; in the latter case on may resort multiple branches, 
or to a parametric two dimensional solution that is found through a Hamiltonian 
approach. Backward propagation is now easily handled, showing that light can follow 
very efficiently a full circle bend waveguide constituted of LLs. Both types of solutions, 
viz. quasi-2D as well as true 2D give the same results in all cases compared. The 
ray tracing results were compared with in house FDTD as well as commercial code 
(COMSOL) and found good agreement. For gaussian beams the LL focusing point 
shifts, an artifact that can be corrected through nonlinearity also in the LL networks. 
The LL networks may find usueful applications in stronly focusing systems [25] . 
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Appendix 

For completion we provide in this Appendix two alternative ways to derive the LL 
ray tracing solution that are based a parametric and a Helmholtz equation approach 
respectively. 

Appendix A. Parametric 2D ray solution 

We use the infinitesimal arc length ds = \J dx 2 + dy 2 in Cartesian coordinates and 
further introduce the parameter t as generalized time; we have ds = \J x' 2 + y' 2 dt with 



a' = da/dt and x = x(t), y = y(t) [H El HU [121 HS1 HZl HH [22]. The Fermat integral 



becomes 
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where n(x, y) the refractive index in cartesian coordinates; the optical Lagrangian 



is 



C(x, y, x', y', t) = n(x, y) a/ x' 2 + y' 2 (A.2) 

We introduce the momenta k x , k y that are conjugate to x, y: 

dC nx' , . , 

k x = TT- = ; (A.3) 
dx' ^/x' 2 + y' 2 V ' 

_ dC _ ny> 

V dy> ^/ x ' 2 + y > 2 1 ' ; 

The Eqns. fl A. 3[) . (I A.4[) consist of an algebraic nonlinear system with solution: 

k 2 x + k 2 y -n(x,y) 2 = (A.5) 
We rewrite Eq. ( 1A.5I) in vector using r = (x, y) and k = (/c^, A;^). Subsequently 

fc 2 - n(r) 2 = (A.6) 
We introduce a new function K(f, k) by 

AT(r,jfe) = ]^{k-n(r) 2 ) = (A.7) 

that has zero total differential, viz. 

dK dK - 

dK(r, k) = — - dr + — ■ dk = A.8 
dr Qk 

We may obtain a Hamiltonian ray tracing system by solving the Hamilton's equations 
for Eq.jEZD [SSlEl] 

dr dK r , A . 

— = — = k A.9 

* = - aF = "2 Vn(rT (A10) 

where V = (Jj, J^) and t is an effective time which is related to the real travel time as 
dt = cdr. Combining finally Eqns. (lA.9p(IA.10p we obtain 

c 2 

r = — Vn(f) 2 (A.ll) 

where in Eq. (lA.llj) we take the derivatives respect to travel time r thus q = dq/dr for 
arbitrary q(r) [21 El EE El EEE1 H5] - Equation (lA.lip is a general equation for arbitrary 
refractive index functions n(r). The explicit solution for LL will be given in the next 
section. 
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Appendix B. Helmholtz wave equation approach 

The stationary states of an monochromatic EM wave are given by the solutions of the 
following Helmholtz equation: 

[\7 2 + (nk ) 2 ]u(x,y)=0 (B.l) 

where u(x, y) is a scalar field, n is the refractive index that generally depends on position, 
k = uj/c = 2ir/\, where to and A are the angular frequency and wavelength of the EM 
wave and c is the velocity of light [2j |3j |T7| [23] . 
Using the well known transformation 

u(x, y) = A(x, y)e l(t>{x ' y) (B.2) 

where A, <fi are real, into the Helmholtz equation ( IB. II) we obtain 

(V0) 2 - (nhY = ^ (B.3) 

V • (A 2 V0) = (B.4) 

The last equation of the above system, Eq.( lB.4l) . express the constancy of the flux 
of the vector A 2 V0 along any tube formed by the field lines of the wave vector defined 
through k = V0; the latter definition turns Eq.( 1B.3l) into 

k 2 - (nk f = ^ (B.5) 

Using r = (x, y) with the operator V = = (J^, ^) and introducing the function 
D(f,k) [23]. 

X7 2 A 



k 2 - (nkn) 2 - 



A (R6) 
we end with a Hamiltonian ray tracing system with the following Hamilton equations: 



df 
It 
dk 
~dt 



dD _ ck 
dk k 
dD 



df 



V 



ck 2 

n 

2 



V 2 A^ 



2k A 

The second term of the Eq. (lB.8l) is the Helmholtz Wave Potential 
W(r) = C ^ A 



(B.7) 
(B.8) 

(B.9) 



2k A 

The potential of Eq.( ]B.9j) keeps the wave behavior in the ray tracing, i.e. the 
diffusion of the beam [23] . In the case when the space variation length L of the beam 
amplitude A(f) satisfies the condition k^L » 1 i.e. Ao << L, the Helmholtz potential 
of Eq. (IB.9j) is approximately zero, thus the Eq.(|R4]) gives the the well known eikonal 
equation, which is the basic equation in geometrical optics approach [7tlT3 |fT7lfT8lfT9ll23] . 

(V0) 2 =(nk ) 2 (B.10) 
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The most important result of this approach is that the rays are not coupled any more 
by the Helmholtz Wave Potential and they propagate independently from one another. 
Finally the equations of motion take the form 
df c r 

— = —k B.ll 

dt k 

d ± = *Vn 2 (B.12) 
dt 2 v ' 

which are described by the Hamiltonian 

H{r,k) = ^P- C ^n\r) (B.13) 
Furthermore the system of equations of motion can be written in a second order 

ode EEjiminidgjdg. 

r 2 

r = — Vn 2 (B.14) 

Using now the LL refractive index of Eq.flT]) we obtain the following equation of 
motion which gives the ray paths inside a LL. 

•■ c 2 

r + -^r = (B.15) 

The solution of Eq. flB.16l) is simple; using boundary conditions r(0) = r~o = (xo, yo) 
and fo = ko = (ko x , ko y ) we obtain: 

ssj-UHMt)?-^ (Bie) 

The solutions ( IB.16j) in Cartesian coordinated describe elliptical orbits, in 
agreement with Luneburg's theory as well as Eq. (1141) . 
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